Journal of Modern Mathematics Frontier Volume 2 Issue 3, September 2013 



www.sjmmf.org 



The Concentration and Dispersion Coefficient 
of the Advection Dispersion Equation 

Tieliang Wang 1 , Kan Zhang* 2 , Xue Feng 3 

^Postdoctoral of Agricultural Engineering, Shenyang Agricultural Uin versify, Shenyang, 110866, China 
2 < 3 College of Sciences, Shenyang Agricultural University, Shenyang 110866, China 
^angzhangpaperfgi^S.com; * 2 kanzhang2004@163.com; 3 xfeng2000@163.com 



Abstract 

In this paper, some recent developments in the analysis of 
the BMO space have been applied to the advection 
dispersion equation. Then the upper and lower bound have 
been obtained for the concentration and dispersion 
coefficient of the advection dispersion equation under the 
conditions where source term and stirring term of advection 
dispersion equation belong to BMO space. 
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Introduction 

Advection dispersion equation is a kind of basic 
equations, which can be used for convection diffusion 
problems. For the passive scalar, complicated behavior 
is often observed even for laminar velocity fields. This 
is the well-known effect of chaotic advection in (Aref 
1984; Ottino 1989). Thus the source and stirring term 
of the advection dispersion equation can be selected to 
be any divergence-free, possibly time-dependent flow 
field. The upper and lower bound of concentration 
and dispersion coefficient play an important role in 
comparing the relative effectiveness of various 
scenarios and analyzing adsorption and desorption of 
some elements in research area, so effort is focused on 
researching the concentration and dispersion 
coefficient of advection dispersion equation in 
(Griecoa 2005; Shawa 2007). It has been confirmed that 
the location of the bounding parabola depends on the 
structure of stirring term and the source term in 
(Plasting 2006). Bounds on the derivative moments of 
the concentration field in (Schumacher 2003) have 
been produced, as well as the bound on mixing 
efficiency for the advection dispersion equation in 
(Thiffeault 2004). 

The classical BMO space was introduced by John- 
Nirenberg in (John 1961; Garnett 1981). A new type of 



BMO, denoted by BMO ' 3 (Q) , has been introduced in 
(Bekolle 1990; Berger 1998) for any bounded domain 
Q of the complex space C" . In this paper, the BMO 
space is defined by the Bergman metric. For 1 < p < oo , 

we write ||| for the norm on L p (B n ,dv a ) and ( v ) a 
for the inner product on L 2 (B n ,dv a ) , where B n is the 
unit ball in C and v a is the Lebesgue volume 
measure on B n . For z&B n , y/ z is the automorphism of 
B n such that y/ z (0) = z and y/ z = , which is 
described in (Rudin 1980). The Bergman space A^{B n ) 
is the space of holomorphic functions which are 
square-integrable with respect to measuring dv a on 

B n . Reproducing kernels K% and normalized 
reproducing kernels k% in A^{B n ) are given by, 
respectively, 

n+a + l 

for z,coe B n . 

For feL(B n ,dv a ) , the Berezin transform of / is 
defined as the function / , that is 

f{z) = \ B J{co)\k a z {co)^dv a {co) 

Let feL(B n ,dv a ) an p>\ , it is said that 
feBMO p (B n ) whenever 

\f\ B MO> =SUp||/°^-7(z)|| <oo 

By the Theorem 5 in (Zhu 1992), the fact is well known 
that BMO p (B n ) is equivalent to BMO p (B n ) in (Zhu 
1992). 
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For z,caeB n/ let 



System Description 



a, ^ l ^ 1 + KH 
)g(z,fl>) = -lQg Z 

2 i-k ? <» 



denote the Bergman metric on B„ . For any zs B n 
and r > , let 

D(z) = {ft)eB n |^(z,tt))<r} 

be the Bergman metric ball with center z and radius 
r . Let \D(z)\ denote the volume of D(z) , i.e. 

\D(z)\ = v a (D(z)). 

For / e l}(B n ,dv a ) , the average of / over D(z) is 
defined by 

1 



f(z)-- 



P(z) 



l D(z) f(ca)dv a (o)) 



Using the properties of Bergman metric, it follows that 



sup 

zeB„ 



1 f I " \P 

} D(z) \f{(o)-f{z)\ dv a (co) <oo 



\D(z) 



if and only if \\f\\ BMO ,, < <*>. 

So the functions in BMO p (B n ) have bounded mean 
oscillation in the Bergman metric and the BMO space 
plays a significant role in operator theory and applied 
science, for example (Duraiswamy 2009; Zhang 2011). 
Since BMO p {B n ) functions are locally mL p {B n ,dv a ), it 
is easy to see that 

^(B n )^BMO p (B n )^L p (B n ,dv a ), p>\ 

BMOl (B n ) c BMO p a (B n ) c BMO x a (B n ), \<p<q 

From now on we will write BMO p instead 
oiBMO p {B n ). 

In this paper, some recent developments in the BMO 
space have been applied to the advection dispersion 
equation. Then the upper and lower bound have been 
obtained for the concentration and dispersion 
coefficient of the advection dispersion equation under 
the conditions where source term and stirring term of 
advection dispersion equation belong to BMO space, 
further the upper and lower bound of the equivalent 
dispersion coefficient are acquired. 

Throughout the paper, the letter c is utilized to denote 
a generic positive constant that can change its value at 
each occurrence. 



The advection dispersion equation for the 
concentration 6{x, t) of a passive scalar is 



d0_ 
dt 



+ D-V0 = kA0 + s$ + sf 



(2.1) 



Where T = top width of the channel; A = wetted cross- 
sectional area; C = area's average dissolved heavy 
metal; Q L = lateral inflow or outflow discharge; 
C L = lateral inflow or outflow dissolved heavy metal 
concentration; Ax = distance between two consecutive 
cross sections which can be either constant or 
variable; H = averaged depth over the cross section; 
U* =local shear velocity; k =dispersion coefficient, 
k > ; u(x,t) = cross-section average velocity; 



D = [7.428 + 1.775(-) a620 (— ) 0572 ]// M (— ) 
H u U» 



5q = — Ajc ; Sf = kCA ; Q, T, H, U* are constants. 

Q 

Without loss of generality, it may be supposed that 
Sq , Sf , 8(x, t) have zero spatial mean at all the times 
and belong to BMO space, which is true under the 
realistic condition. Then research domain is divided in 
which finite unit ball B n is taken into consideration 
and at the same time, it is assumed that both the 
source term and stirring term act on a comparable 
scale \ <l< L ,where L is length of side of research 

area and — is an integer. The research time is 
supposed to be — on the unit ball, where v is an 

T 

appropriate time scale characterizing the source and 
stirring term and 0<f<l . For 

< x <l, S$ (-,-), Sf (-,-) and «(-,-) are the source 
It It It 

term and stirring term on the unit ball. In order to 

keep track of the effects of the amplitudes of the 

source term and stirring term on the upper and lower 

bound for the concentration ||#|l , and dispersion 

II IIBMO; r 

coefficient k , supposing that Sq + s f\ = $ an d 

II Wbmo,, 



WBMOf, 



U. 



Suppose the advection dispersion operator defined by 

L = — + DV-M 

dt 

then the adjoint of the advection dispersion operator is 
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L = — -D-V-kA 
dt 

Next, the upper and lower bound of 
IHIbmo 2 W1 ll be obtained for given 5 and U . 

Since 



S ~ IK + S < \bMOI ~ W L \bM01 P\\bM01 ~ \\ L \\bM01 ^BMOl 



then we have 



S< sup [ 

IMLo2=i 



dy/ 
~dT 



-D-Sly/ -kAy/ 



bmo;, 



PWbmoI (Z2) 



Where y/ is an arbitrary smooth noralized spatially 
periodic function. 

By the application of the Minkowski inequality and 
Cauchy inequality to formula (2.2), it can be observed 
that 



S < [ sup 



\\BMO a =1 



dy/ 



dt 



+ k sup \Ay/\ 



BMOi 



BMOi 



\\BMO a =\ 



+ D 



BMOt 



BMOi J H \\BMO;, 



WBMOa =1 



< [ sup 



dy/ 



dt 



-k sup \\Ay/\ 



7.428// m 



+ - 



BMO L a 
2 

BMO: 



BMOi 



U. 



— sup \wy/\ 



BMOi 



(2- 



, 0.38^0.6211 ||l-428 

1.775// T m 2 

+ ^ sup v^l .M 



(U.) 



0.428 



n. 



WBMOi J ll llBMOZ 



By the formula (2.3), it is clear that 



WBMOz 



cU 



(2.5) 



By the formula (2.5), it is clear that we have the 
freedom to choose y/ to optimize the lower bound of 
the concentration for the particular source and stirring 
term. From formula (2.5), the fact is learnt that 
increasing U at fixed S should decrease the lower 
bound of concentration. While increasing S at fixed 
U should augment the lower bound. Conversely, the 
increasing of concentration does not imply the 
increment of source amplitude or stirring amplitude. 



However, the concentration ||#|| 



with large 



amplitude almost produces the large source amplitude 
S unless U is decreased sufficiently. It will be proved 
that the lower bound of the concentration is 
SI 

proportional to — , so that a source with large 

fluctuations necessarily produces a poorly mixed state 
unless U is increased sufficiently. 

By formula (2.4), (2.5) and substituting the scaled 



variables T = — and y ■■ 

T 



WBMOi 



I 



, it is easy to prove that 



SI 



U[S r sup 



dy/ 
dT 



BMOi 



sup \\A v y/\\ . +c, 
,ii r II >r WbmoI 



(2.6) 



Where 



:C HI fibril BMO;; 



= cU\\9\ 



Where 



BMOi 



sup 



dy/ 



dt 



BMO- 



C = - 



+ k sup \\Ay/\ 



BMOi 



WBMOi 



7.428// m 



BMO; 



u. 



— sup \Wy/\ 



BMOi 



1 -7-7<rr/ 0.38^0.620 II ||0.428 

II \\BMO„ II,-, | 

+ TTT^-r, — SUp \\Vy/\ 



0.428 



ML 



BMOi 



3) 



(2.4) 



7.428// | 
q=— sup \\S7y/\ 



BMOi 



1 775// °' 38 7"°- 620 

~ ,0.428,1 ||0.572 , SU P ll V ^IL,/ ( , 



(2.7) 



wu\ 



bmo; 



' BMO a -1 



Here P=—is the Peclet number and S r = — is the 

e k r Ut 

Strouhal number. 

If it possible to choose y/ to be time-independent and 
still satisfy formula (2.6) and (2.7), then we get the 
lower bound of concentration 



II WbmoI 



SI 



U[P e sup \\A x .y/\\ 2 +c,] 
c i, ,i II > T WbmoI 



(2.8) 



WBMOa =1 
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By the formula (2.6) and the assumption that y/ is 
time-dependent, then the lower of concentration is 
affected by Peclet number and Strouhal number, 
holding the other parameter constant. By formula (2.8) 
and the assumption that y/ is time-independent, then 
the lower of concentration is only affected by Peclet 
number, holding the other parameter constant. As 
known, the concentration and dispersion coefficient 
can be affected by several environmental and 
biochemical factors such as temperature, PH, EC, 
salinity, etc (Nassehi 1993; Kashefipour 2002; 
Balmforth 2003; Roshanfekr 2008). In fact, formula (2.6) 
and (2.8) also imply that there must have lots of 
complex factors influencing the concentration and 
transfer theory. 

By the formula (2.4) and (2.5), we get the lower bound 
of dispersion coefficient k . 



k> 



sup 

'IL |, 



8y/ 



BMO; 



fWbmoi „ ,? up IMU>; 

WhuoU 



7.428// u 



, SU P \M\bmoI U -\H 



sup ||V(i/| 



sup A;// 



bmo;, 



bmo;, 



(2.9) 



1 ,|| ||2 

x-y = — ( \\x + y 
4 11 



BMO„ 



for x,y e BMO 2 , it is clear that 



4 +sf - — + kA0 
dt 



BMO-;, 



= \\D-V0\ 



BMO' 



^P\Lol + \^<MOl-\\ D -^ d tBMOl) 
^\i D tBMOl +d P\Lol-\ D -^<MOl) 



where d is poincar e s number. 

By the formula (2.10), it is easy to prove 



ii ii 2 ii ii 2 

\\ D \\bM01 +d \\ \\BMOl 



<2 



sH+s?- — + kA0 



dt 



+ D-ve 



MIO; 



BMOt 



<2S + 2 



80 



dt 



BMO;, 



+ ^MbmoI + \\ D \\bmoI 



+ A\ d \\bmo^ 9 \\bmoI + W 9 \bmoI 



By formula (2.11), it is obtained that 



(2.10) 



(2.11) 



1 77<IW - 38 T a620 IUI L428 

1.775// T \u\\ mol 



sup \Wy/\ 



7P IMUog^-^ikUa. 



bmo;, 



If y/ is time-independent and satisfies the formula (2.9), 
then we have 



k> 







bmo;. 



sup Ay/ - 



bmo;. 



7.428/7 m 



BMOf x =1 

2 



\bmo: 



sup Ay/ - 



bmo; 



^77 „ , su p Wn BM oi 



n 



dm 



<2S + 2 



<2S + 2 



cO 



■^M\ B M0l + A D \ B M0lW d \\ B M^ 



\\V0\ 



BMOl 



at 



7.428// Ml 



+ n A0 \\ B MOl + \\ V0 \\ B MOl + ^- 



U„ 



<25 + 2 



+ ^M B M0l + W d \\ B M0l +( -- 



14.856//|w||' 
U, 



3.55r"-//0.3 8[/t -o.4 2 s |H| ^ ||V , L 



(2.12) 



1.775// a38 r - 620 || M ||^ 2 



sup V y/\ 



bmo;, 



Let 



30 



dt 



BMO; 



bmo; 



Next, we will obtain the upper bound for ||#| 



bmo; 



Using poincare's inequality, Frechet-Jordan-Von 
Neumann theorem and the fact that 



By formula (2.12), we obtain 

\\0\\ BMO 2 < c^2S + 2s x + 2ks 2 + (c 3 U - 512 + c 4 )U 1A2S s 3 + £ 



(2.13) 



88 



Journal of Modern Mathematics Frontier Volume 2 Issue 3, September 2013 



www.sjmmf.org 



H 



^0.38^0.620 



c 3 = 14.856 — ,c 4 =3.55 

U. (U,) 0A2S 

If 6 is time-independent and still satisfies formula 
(2.13), then we acquire 

\\e\\ BM0l < c^2S + 2ks 2 + (c 3 U - 512 + c 4 )U lA2 *s 3 + si (2.14) 



^0.38^,0.620 



c 3 = 14.856 — ,c 4 =3.55 , ... 

£/. (i/,) 0428 



By the formula (2.9) and (2.16), we obtain the upper 
and lower bound of the dispersion coefficient k . 



S 



sup 

IML 2=i 



d\j/ 



BMOi 



sup I A ^| 



BMO„ 



WBMOa =1 



7.428// U 



sup 



BMO„ 



By the formula (2.5) and (2.13), we get the upper and 
lower bound for the concentration lldll 



WBMOi 



cU 



WBMOi 



< cyjlS + 2e 1 + 2ke 2 + (c 3 U°- 512 + c 4 )U lAn s 3 + s 



(2.15) 



If 9 is time-independent and still satisfies formula (2.8) 
and (2.14), then we obtain 



SI 



U[P e l sup A V J , +c] 



< \\9\ 



BMOi 



<c^2S + 2ke 2 +{c 3 U Q - 512 +c 4 )U 1A2S e 3 +e 



From the formula (2.15), the consequence of the bound 
of upper and lower bound for ||#|| BM0 2 is that large 

IHIbmo 2 must eventually imply large S , holding the 

other parameters constants. However, the large S 
does not necessarily imply large ||#| BM0 2 • 

Using the same method, we obtain the upper bound of 
dispersion coefficient k . 



k<[S + 



80 



dt 



+ 1.5x 



BMOi 



(7.428— U 2 +1.7757/0.38^0.620^-0.428^1.428 (2 16) 



+ W d \\BM 9 H^lsMOl 



If 6 is time-independent and satisfies the formula 
(2.16), then we have 

k<[S+1.5x 

(7.428— U 2 +1.775// a38 7; a620 [/,; - 428 [/ 1 - 428 +||V^|L _ 2 ) 2 ] 



Ad 



<k<[S + 



(7.428 — U 2 + 1.775//0.38 :r 0.620 f/ -0.428 f/ 1.428 

u. 



L77 5 H"«r^>|»g, 
5(9 



(2.17) 



dt 



+ 1.5x 



BMO, f 



+ wn B M iyH^\\BMoi 



By the formula (2.17), the increment of U should 
decrease the lower bound of dispersion coefficient and 
increase the upper bound of dispersion coefficient, 
holding the other parameters constants. While 
increasing S should augment the upper and lower 
bound of dispersion coefficient, holding the other 
parameters constants. Conversely, the increment of the 
upper and lower bound of dispersion coefficient does 
not imply the increment of source amplitude or the 
decrement of stirring amplitude. However, the upper 
and lower bound of dispersion coefficient with large 
amplitude almost produce the large source amplitude 
S , unless U is decreased sufficiently. 

As a physically meaningful measure of mixing 
efficiency, the quivalent dispersion coefficient k eq has 

been investigated in (Young 1999; Pope 2000). Let 
6q be the solution of the advection dispersion equation 
with the same source term but no stirring term, then 
we define 



0n 



\BMO- 



\o\ 



(2.18) 



BMOf, 



By the formula (2.17) and (2.18), it is easy to obtain the 
upper and lower bound of the equivalent dispersion 

coefficient k 2 q . 



BMOf, 
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sup 



dy/ 



BMOi 



PWbmoi, ,f up IMUoj 

HI, 



sup 



BMO:, 



\\BM0a =1 



7.428// U 



mo,: 



1.775// T \\u\\ BM0l 



0.428 



sup \Wy/\ 



BMO;, 



>\\BW%_ 2 
,,2 - K eq 

WbmoI 



(2.19) 



II IIbmo„ 



01 



+ 1.5x 



BMO„ 



(7.428 — (7 2 + 1 . 775// 0.38 r 0.62 0(/ -0.42 8(/ 1.428 



+ HUp 2 ] 2 HMI 2 2 



The worst lower bound for the mixing efficiency 
would be achieved by injecting scalar variance at scale 
/ while stirring to keep the dominant scale of the 
concentration fluctuation field as L. 

Conclusions 

By the same method, it is easy to obtain the upper and 
lower bound for the other parameter of the advection 
dispersion equation and at the same time we may find 
the relation among the parameters. It is clear that our 
results are established on the unit ball B n , however, 

our results are obviously correct for the whole 
research area with any bounded domain Q . It is 
assumed that 9 i is the concentration in the ith unit ball 
and kj is the dispersion coefficient in the ith unit ball, 

so 



9- 



l + 2 +--- + n 



is the concentration of the whole research area and 
k l +lc 2 +--- + k„ 



k = 



is the concentration of the whole research area. Using 
the same method, the upper and lower bound are 
obtained for the concentration and dispersion 
coefficient for the whole research area, as well as the 



upper and lower bound of equivalent dispersion 
coefficient for the whole research area. 

The theoretical results obtained in this paper will be 
applied to the realistic examples, for example, to the 
research for the transfer of dissolved heavy metal in 
river. According to the formula (2.15), the following 
parameters are required: top width of the channel; 
wetted cross-sectional area; the average concentration 
of dissolved heavy metal for the research area; average 
depth over the cross section; cross section average 
velocity; local shear velocity and the selection and 
calculation of iff . It is clear that most of these 
parameters are easy to be obtained by adopting 
general method, so new way to research the 
concentration for the dissolved heavy metal is 
available. 

The investigation of these works is left for further 
research. 

REFERENCES 

Aref, H. "Stirring by chaotic advecton." Journal of 
Fluid Mechanics 143(1984): 1-21. Balmforth, N. J., and 
Young, W. R. "Diffusion-limited scalar cascades." 
Journal of Fluid Mechanics 482(2003): 91-100. 

Bekolle, D., Berger, C. A., Coburn, L. A., and Zhu, K. H. 
"BMO in the Bergman metric on bounded symmetric 
domains." Journal of Functional Analysis 93(1990): 310- 
350. 

Berger, C. A., Coburn, L. A., and Zhu, K. H. "Function 
theory on Cartan domains and the Berezin-Toeplitz 
symbol calculus." American Journal of Mathematics 110 
(1998): 921-953. 

Duraiswamy, K., Manjula, D., and Maheswari, N. 
"Advanced Approach in Sensitive Rule Hiding." 
Modern Applied Science 3(2009): 98-105. 

Garnett, J. Bounded analytic functions. Academic Press, 
New York, 1981. 

Griecoa, L., Tremblayb, L. B., and Zambianchia, E. "A hybrid 
approach to transport processes in the Gulf of Naples: an 
application to phytoplankton and zooplankton 
population dynamics." Continental Shelf Research 
25(2005): 711-728. 

John, F., and Nirenberg, L. "On functions of bounded mean 
oscillation." Communications on Pure and Applied 
Mathematics 14(1961): 415-426. 



90 



Journal of Modern Mathematics Frontier Volume 2 Issue 3, September 2013 



www.sjmmf.org 



Kashefipour, S. M., and Falconer, R. A. "Longitudinal 

dispersion coefficient in natural channels." Water 

Research 40(2002): 1596-1608. 
Nassehi, V., and Bikangaga, J. H. "A mathematical model for 

hydrodynamics and pollutant transport in long and 

narrow tidal rivers." Applied mathematical modelling 

17(1993): 415^22. 
Ottino, J. M. The kinematics of Mixing: Stretching, Chaos, 

and Transprot. Cambridge University press, 1989. 
Plasting, S. C, and Young, W. R. "A bound on scalar 

variance for the advection-diffusion equation." Journal 

of Fluid Mechanics 552(2006): 289-298. 
Pope, S. B. Turbulent flows. Cambridge University Press, 

2000. 

Roshanfekr, A., Kashefipour, S. M., and Jafarzadeh, N. 
"Numerical modelling of heavy metals for riverine 
systems using a new approach to the source term in the 
ADE." Journal of Hydroinformatics 10(2008): 245-25. 

Rudin, W. Function Theory in the Unit Ball of C" , Springer- 
Verlage, New-York, 1980. 

Schumacher, J., Sreenivasan, K. R., and Yeung, P. K. 
"Schmidt number dependence of derivative moments for 



quasi statie sraining motions." Journalof 
Fluid Mechanics 479(2003): 221-230. 
Shawa, T. A., Thiffeaultb, J. L., and Doering, C. R. "Stirring 
up trouble: Multi-scale mixing measures for steady scalar 
sources." Physica D: Nonlinear Phenomena 231(2007): 
143-164. 

Thiffeault, J. L., Doering, C. R., and Gibbon, J. D. "A bound 
on mixing efficiency for the advection-diffusion 
equation." Journal of Fluid Mechanics 521(2004): 105- 
114. 

Young, W. R. Stirring and Mixing: Proc. 1999 Summer 
Program in Geophysical Fluid Dynamics. Woods Hole 
oceanographic institution, Woods Hole, MA, USA, 1999. 

Zhang, K, Liu, C. M., and Lu, Y. F. "Toeplitz operators with 
BMO symbols on the weighted Bergman space of the 
unit ball." Acta Mathematica Sinica (English Series) 
27(2011): 2129-2142. 

Zhu, K. H. "BMO and Hankel operators on Bergman 
spaces." Pacific Journal of Mathematics 155(1992): 377- 
395. 



91 



